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Abstract. New ages are computed for the stars from the 
Edvardsson et al. (1993) data set. The revised values are 
systematically larger toward older ages (t>4 Gyr), while 
they are slightly lower for t < 4 Gyr. A similar, but con- 
siderably smaller trend is present when the ages are com- 
puted with the distances based on Hipparcos parallaxes. 
The resulting age-metallicity relation has a small, but dis- 
tinct slope of ~0.07 dex/Gyr. 
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1. Introduction 

The Edvardsson et al. (1993, hereafter Edv93ea) sample 
of F- & G-stars show a slowly, increasing gradient with 
decreasing age in the age-metallicity relation (AMR, see 
Edv93ea Fig. 14a). There is mainly a large scatter in 
metallicity at a particular age or vice versa. This scatter 
could be real or it resides cither partly in the metallicity 
or in the age determined for each star. 

We focus the present analysis on the age. A bias in the 
ages could originate from the isochrones or the method 
adopted for the age determination. Edv93ea used the 
Vandenberg (1985, hereafter Vdb85) isochrones to deter- 
mine the age in the (log T e g, AM V ) diagram. We used the 
isochrones from Bertelli et al. (1994, hereafter Bert94ea) 
to determine the age in the (log Toff, My) diagram, see 
Sect. 2.2.2 for details. Moreover, the new distances from 
Hipparcos parallaxes (ESA 1997) further justify a re- 
examination of the stellar ages of the Edv93ea sample. 

The aim of this paper is to determine new ages for the 
stars from the Edv93ea sample. Section 2 gives a brief de- 



* Tables 3-8 are available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79,128,5) or 
WWW at URL http://cdsweb.u-strasbg.fr/Abstract.html 



scription of the data set, an outline of the method used to 
determine the ages of the stars and the results. In Sect. 3 
we present the newly obtained AMR and discuss the lim- 
itations of the analysis. 

2. Analysis 

2.1. Data 

The sample comprises mainly a volume limited (d < 50 pc) 
set of F- & G-type stars for which the distances are de- 
termined (Edv93ea - photometric distances; ESA 1997 - 
trigonometric parallaxes) and the effective temperatures 
& metallicities are determined from spectra (see Edv93ea 
for details). The V-band photometry of these stars are 
taken from Gr0nbech & Olsen (1976), Olsen (1983, 1993), 
Perry, Olsen & Crawford (1987), Schuster & Nissen (1988), 
and Hoffleitfe Warren (1991). 

2.2. Method 

2.2.1. The isochrones 

We used the Bert94ea isochrones for the determination of 
stellar ages from the Edv93ea sample. The initial chemical 
composition of the isochrones parameterized with (Y, Z) 
obeys the relation AY/AZ = 2.5 (Pagel 1989). The grid of 
isochrones is obtained from stellar models computed with 
the radiative opacities from OPAL (Iglesias et al. 1992, 
Rodgers&Iglesias 1992). They have (Y,Z): (0.23,0.0004), 
(0.24,0.004), (0.25,0.008), (0.28,0.02), and (0.352,0.05). 
One isochrone grid with (Y, Z) — (0.23, 0.001) is based 
on models calculated with the radiative opacities from 
Huebner et al. (1977). This complete set of isochrones 
of fixed metallicity is interpolated to obtain isochrones 
with intermediate metallicities. The metallicity range, 
smoothly covered by the isochrones, thus spans from 
Z = 0.0004-0.05. 

This enables us to generate isochrones with metallici- 
ties comparable to those used by Edv93ea. There are how- 
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Table 1. Sensitivity of the age determination through varia- 
tion of the input parameters (effective temperature, metallicity, 
distance). Bias corrections (Lutz&Kelker 1973, see Sect. 2.2.5 
for additional details) were not applied. For each parameter 
with a good fit to the data the age difference 8 is calcu- 
lated with respect to the mean age obtained for the star, i.e. 
8 — log(age„) — log(agc). The table gives the mean age dif- 
ference 5 for each parameter together with its variance a. The 
first set demonstrates the sensitivity of the parameters with 
the uncertainties as discussed in Sect. 2.2.3. In the second and 
third set of calculations we assume respectively a 10% and 5% 
uncertainty in the distance. In the last set we assume also a 
smaller uncertainty for T e ff and [Fe/H] for a 10% uncertainty 
in the distance 
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over marked differences between the Vdb85 and Bcrt94ca 
isochrones. They are respectively: 

— the helium fraction, fixed Y = 0.25 versus variable 
through AY/AZ = 2.5; 

— the radiative opacities, LAOL (Huebner et al. 1977) 
versus OPAL (Iglesias ct al. 1992, Rodgers & Iglesias 
1992); 

— initial metal abundance mix of elements heavier than 
helium, Vandenberg (1983) versus Grevesse (1991); 

— convective overshoot, none versus included. 

In addition, different tables have been used by Vdb85 and 
Bcrt94ca to convert the effective temperature and lumi- 



Table 2. Similar to Table 1, but now with the distances de- 
rived from the Hipparcos parallaxes (ESA 1997). 
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nosity into magnitude and colours. The Vdb85 isochrones 
were computed for metal mass fractions of Z = 0.0169, 
0.01, 0.006, 0.003, and 0.0017. Note that for the conver- 
sion of [Fc/H] or [Mc/H] to Z the reference solar metal- 
licity is Z (Vdb85) = 0.0169 and Z (Bert94ea) =0.020. 
In addition the Vdb85 isochrones ought to be shifted 
by (51ogT ef j = -0.009 to satisfy the solar constraint (see 
Edv93ca pi 19 for details). Figure 1 shows a set of 
[Fe/H] =0.0 isochrones from Vdb85 and Bert94ea. Note 
the different appearance of the 10 Gyr isochronc from 
Bcrt94ea with respect to isochrones of younger age due 
to the presence of convective overshoot in the latter. 

2.2.2. Fitting (log T off , My) 

For each star in the sample the metallicity is known. 
Each star can be located in a Hertzsprung- Russell diagram 
(HRD) through its effective temperature and absolute vi- 
sual magnitude. We used the interpolated isochrones from 
Bcrt94ca to determine the age of the star, where the stel- 
lar mass specifies the position along the selected isochrone. 
Powell's method (Press et al. 1986) is used to obtain the 
stellar age and mass. We minimize the difference in the 
(logT e ff, My) diagram between the observed star and the 
interpolated isochrones from Bert94ea (see for example 
Fig. 1). For this purpose we defined a weighted, reduced 
chi-squared 



y (My,„ - My.jJ +16(^l0gT efi> - logTeff,^ 

The subscripts o and i refer respectively to the ob- 
served and the synthetic, isochrone quantities. We used 
the Bert94ea relation — log Z/Z Q = 0.977[Fe/H] — to 
convert [Fe/H] to Z. 

Note, that our fitting method is not substantially 
different from the determination of the age in the 
(logT e ff, AMy) diagram used by Edv93ea. Since AMy Q 
= My QiZams - My and AMy ;i = M v ^ zams - M Vii , then 
minimizing the difference between the observed and the 
isochrone magnitude implies: SMy — AMy,o — AMy,i- 
With Mv,o,zams = Mv,i, Z ams this equals to the quantity 
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Vandenberg (1985) 
Bertelli ct al. (1994) 



3.85 




3.8 3.75 
log T eff 



Fig. 1. Comparison for [Fe/H]=0.0 between the Vdb85, 
LAOL based isochrones with the OPAL based isochrones from 
Bert94ea. The isochrones ages are 2, 4, 6, 8 and 10 Gyr 
for the LAOL based set and 2, 4, 6.3, 8 and 10 Gyr for 
the OPAL based set. The Vdb85 isochrones are shifted by 
<51ogT c ff =-0.009 to satisfy the solar constraint (see Edv93ea 
for details). In addition, the stars (Table 4) with metallicities 
— 0.05< [Fe/H] <0.05 have been inserted. The symbols are the 
same as those in the following figures 
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Fig. 2. A comparison between the ages from Tables 3&4 with 
those determined by Edv93ea, due to different isochrone sets. 
The open circles denote the MS stars from Table 3 and the 
asterisks denote the SGB stars from Table 4. The big symbols 
are for stars with at least 6 out of the maximum 7 good age 
obtained from the analysis (see Sect. 2.2.3 for details) and for 
which the uncertainty in log(age) is smaller than 0.05. The 
small symbols are used in all remaining cases and refer in gen- 
eral to stars with less reliable ages 



minimized in the reduced chi-squared denned above. If on 
the other hand My jD>zams ^ My^ zams then an unneces- 
sary bias would have been added in the results. In that 
case minimization in the (logT c g, My) plane should be 
preferred over minimization in the (logT e ff, AMy) plane. 



for each photometric measurement. Then an average was 
computed, weighted over the photometric errors given by 
the various authors or adopted by us (0" l 01 for the pho- 
tometry from the bright star catalogue and 0" l 05 if the 
photometric errors were not specified at all). 



2.2.3. Uncertainties in T off , [Fe/H] & distance 

The uncertainties in the effective temperature, metallicity 
and the distance will give different results for the fitted 
mass and age along the isochrone. For each star we re- 
peated the calculations with slightly different values for 
the effective temperature (-50/+100 K; see Edv93ea), the 
metallicity (±0.07 in [Fe/H]; adopted) and the distance 
(±15% for distances given by Edv93ea; in Tables 7&8 the 
uncertainties are considered separately for each star with 
an Hipparcos parallax). In total 7 different values for the 
mass and the age of each star are obtained. A consequence 
of the automatic fitting procedure is that it always pro- 
vides an answer, which sometimes is not a good fit to the 
actual data. The latter were not considered in the calcu- 
lation of the unweighted mean and the uncertainty in the 
age and mass of each star in the sample. 

When various measures for a particular star of the 
V-band photometry were found in different databases 
the complete procedure as described above was applied 



2.2.4. Isochrone population 

The method thus far considers only the formal solution, 
automatically generated by the procedure outlined above. 
However, ambiguities in the age determinations are not 
taken into account when the stars are located near the end 
of the core-H burning phase. In the vicinity of the termi- 
nation point (i.e. the hook feature) of young isochrones, 
it is possible that the star is either at the end of the main 
sequence phase or beyond the termination point, each im- 
plying a different value for the age of the star. 
In the mass range from 1.0-1.7 Mq the evolution along 
the SGB (Sub-Giant Branch) is roughly 4-40 times faster 
than the evolution from the ZAMS (Zero Age Main Se- 
quence) up to the termination point. Evolved stars must 
therefore be inspected individually. If the resulting age 
results in a position on the isochrone beyond the termi- 
nation point, then an alternative age corresponding to a 
point below the termination point is assigned to the star 
(see last column Table 4). 
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Fig. 3. The AMR for all the stars in Tables 3&4. See Fig. 2 
for the symbol description, but the small symbols are used 
here for the remaining cases with n > 2 and an uncertainty in 
log(age)<0.3 



2.2.5. The Lutz-Kelker correction 

Lutz&Kclkcr (1973) demonstrated that a bias is present 
among stars with their absolute magnitudes obtained from 
trigonometric parallax. The bias is not confined to volume 
limited samples as discussed by Trumpler & Weaver (1953) 
and it depends only on the ratio (ov/V), i.e. the ratio of 
the standard error of the parallax over the parallax. Due 
to observational errors in the parallax one expects to find 
statistically more stars from a larger distance scattering 
into a smaller volume than vice versa. The consequence of 
this geometrical effect is that the stars in a sample are on 
average brighter and younger than they appear to be. 

Thus far we have considered the formal uncertainties 
in the method to determine the age of the stars. Edv93ea 
compared their photometrically derived distances with 
those obtained from ground based parallaxes. They found 
an excellent agreement and estimated that the uncertainty 
in the distances is about 15%. Their a posteriori informa- 
tion about the compatibility with parallax data implies 
that we ought to apply a Lutz-Kelker correction in our 
analysis. This correction is in the generalized case only 
invoked by the relative uncertainties in the distance scale. 
Note however, that Edv93ea (pl21 and references cited 
therein) do not have to apply a Lutz-Kelker correction, 
because the distances that we adopted from their paper 
were neither required or used in their AMy-method. 

The Lutz-Kelker correction is applied to the absolute 
magnitudes determined from distances for the stars from 
Edv93ea. We expect, within the volume considered, a ho- 
mogeneous distribution for the F- & G-type stars and we 
used Hanson's equation 31 (Hanson 1979) with index n = 4 
to compute the value of the correction. The stars become 
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Fig. 4. (a) The AMR for the stars in Tables 3&4 with a disper- 
sion in log(age) <0.05 (see Fig. 2 for the symbol description), 
(b) the AMR for the same stars, but now with the 'visually' 
assigned ages for the SGB stars (filled triangles, last column 
Table 4), and (c) the AMR for the same stars, but now with 
the ages from Edv93ea. 

about 0™3 brighter with cr w /7r~0.15. The correction leads 
to negligible differences in the results with distances from 
the Hipparcos parallax, because for the majority of the 
stars ov /ir <0.03. 

2.3. Results 

The results are presented in three stages. First, we show 
the differences due to a change of isochrones: Vdb85 versus 
Bert94ea. However, no corrections were made for a bias in 
the absolute magnitudes. This is done in the second stage 
of the analysis. In the third and final stage we based the 
analysis on the Hipparcos parallax. 

Table 1 shows the sensitivity of the age through indi- 
vidual variation of each of the input parameters in com- 
parison with the mean age obtained for each star. The first 
set of calculations show the sensitivity of the results for 
the uncertainties, as described in Sect. 2.2.3, adopted for 
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Fig. 5. A comparison of the ages computed with the distances 
of the stars based on Hipparcos parallaxes (Tables 7& 8). Panel 
(a) shows the comparison with the ages from Edv93ea. The 
differences are due to the isochrones and the fitting method. 
In panel (b) we make a comparison with the ages from Ta- 
bles 5&6. The differences are only due the distances adopted 
for the stars (photometric distances from Edv93ea with the 
Lutz-Kelker correction versus distances from trigonometric 
parallax) . Symbols are the same as described in Fig. 2 



the input parameters. In the second & third set of calcu- 
lations a smaller uncertainty in the distance was adopted 
and in the fourth set we further assumed smaller uncer- 
tainties in the effective temperature and metallicity. 
The average results, from the procedure outlined in pre- 
ceding sections, are given in the Tables 3 & 4. They list 
respectively the stars which are on or very near to the MS 
(Main Sequence) and those which are on the SGB. The 
Lutz-Kelker correction was not applied in order to discuss 
the effects from different sets of isochrones. Tables 5 & 6 
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Fig. 6. The AMR similar to Fig. 3, but now for all the stars 
in Tables 7&8 with their distances based on Hipparcos paral- 
laxes (ESA 1997). See Fig. 2 for the symbol description. Small 
symbols are used here for the remaining cases with n > 2 and 
an uncertainty in log(age) < 0.3 



contain the ages when the correction is applied. 
Table 2 shows, similar to Table 1, the sensitivity of the 
age through variation of the individual parameters, but 
now based on Hipparcos distances instead of the distances 
given by Edv93ea. The average results for each star (MS 
and SGB) are given in the Tables 7&8. 

In the following set of figures we first compare the 
difference due to a different choice of isochrones, i.e. 
Figs. 2-4. Figs. 5ab&6 show the relations when the 
trigonometric parallaxes from Hipparcos are adopted. 
In Fig. 2 we compare the ages for the stars computed with 
the Bert94ea isochrones with the ages given by Edv93ea. 
The comparison shows that towards older age the re- 
vised values are systematically larger than Edv93ea, see 
Sect. 3.1 for a detailed discussion about this difference. 
Note that we refer to the stars with reliable ages, i.e. the 
big symbols, and discarded the old MS stars (small open 
circles) from the comparison. 

Figure 3 shows the age-metallicity for all the stars listed 
in Tables 3 & 4. Figure 4a shows the AMR for the SGB 
stars with an average age obtained from at least six out 
of the maximum of seven good age determinations with 
an uncertainty in log(age) < 0.05. Figures 4b & 4c show for 
comparison the AMR for the same stars, but now with re- 
spectively the 'visually' assigned ages for the SGB stars 
(see also Sect. 2.2.4) and the ages from Edv94ea. 
In Figs. 5a & 5b we compare the ages (Tables 7& 8), based 
on absolute magnitudes obtained from Hipparcos paral- 
laxes, with those obtained by Edv93ea and the ages given 
in Tables 5&6. Figure 6 is similar to Fig. 3, but now the 
ages are based on the distances from the Hipparcos paral- 
lax instead of the distances given by Edv93ea. 
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Figure 7 a displays the relation obtained from distances de- 
rived from the Hipparcos parallax, while Fig. 7b shows for 
comparison the same stars, but now with the ages from 
Edv93ea. 

3. Discussion 

3.1. Revised ages 
3.1.1. Isochroncs 

In Fig. 2 a comparison is made between the ages obtained 
by Edv93ea and those determined here with the Bert94ea 
isochrones. The general trend is that the revised ages get 
systematically larger towards older ages. On the other 
hand, the revised ages of a few MS stars are found to be 
considerably younger, while a few other young MS stars 
are now considerably older. 

The general trend should mainly originate from the dif- 
ferences between the Vdb85 and Bert94ea isochrones. Dif- 
ferences in the opacities, initial abundance mix, and the 
different conversion tables have barely influence on the age 
obtained. 

• Convective overshoot is only present in stars with 
M>1 M Q or t^lO Gyr (see Chiosi et al. 1992, Bert94ea, 
and references cited in those papers). The effect of convec- 
tive overshoot at the MS is that the stars have larger cores. 
This results in a higher luminosity and a prolonged age. It 
should lead to systematically older ages for the stars with 
respect to ages determined from the Vdb85 isochrones. 

• The differences in the helium mass fraction has the 
following effect on the age determination: 

— for solar metallicity isochrones the helium mass frac- 
tion Y is lower in the Vdb85 isochrones; the effects 
of difference in Y was checked with the isochrones 
present in our data base with (Z — 0.05,F = 0.352) and 
(Z = 0.05,y = 0.4, unpublished); a lower Y leads to a 
lower T c ff due to an increased opacity from a higher 
hydrogen mass fraction X, while the turn-off luminos- 
ity remains comparable; as a result the MS turn-off 
ages determined with isochrones with a lower Y will 
be younger; from this argument one expects older ages 
from the Bcrt94ea isochrones; 

— for metal-poor isochrones Y is higher in the Vdb85 
isochrones; this leads to the reverse effect, i.e. younger 
ages from the Bert94ea isochrones. 

The general trends outlined above are at young age not 
in agreement with Fig. 2, because we expect older ages 
from the Bert94ea isochrones. Instead we find a compara- 
ble or even a slightly younger age. A similar discrepancy 
is found for old stars with relatively low metallicities. We 
further note, that the Vdb85 isochrones with [Fe/H] =0.0 
were not as expected from the difference in Y located 
at lower, but higher effective temperature. In addition, 
the shift of 5 log T cff = -0.009 applied by Edv93ea to the 
Vdb85 isochrones makes it even more difficult to deter- 
mine the origin of the differences between the Edv93ea 



and our revised ages. 

There are too many differences (e.g. opacities, convec- 
tive overshoot, normalization, conversion tables ...) be- 
tween the Bert94ea and the Vdb85 isochrones that it is 
not only hard to compare them, but also to understand 
which of the differences or a combination of them is/are 
the dominant factor (s). All the factors mentioned above 
hamper the analysis on the origin of the differences be- 
tween the Edv93ea and the revised ages. Although we 
are able to outline general differences between the re- 
normalized Vdb85 and Bert94ea isochrones, we fail to 
prove with plain arguments the exact differences between 
two isochrone sets and hence the differences between the 
ages obtained by Edv93ea and this work. To understand 
the difference between these isochrones requires a deep 
and detailed analysis of isochrones with the same Y, Z. 
Unfortunately they are not available to us, while their 
computation is beyond the scope of this paper. 

3.1.2. Hipparcos parallaxes 

Figure 5a shows the effects due to differences in the 
isochrones and the fitting method. The reliable ages of the 
stars based on the Hipparcos parallax (the big symbols) 
get ~ 2 Gyr older for stars in the age spanning 9-16 Gyr, 
~ 1 Gyr older for stars in the age spanning 4-9 Gyr, and 
a negligible difference for stars younger than 4 Gyr. Figure 
5b displays the effects due to differences in the distance: 
photometric distances from Edv93ea and the derived ab- 
solute magnitudes corrected for the Lutz-Kelker bias ver- 
sus the distances derived from the Hipparcos parallax. It 
indicates that, besides some increased scatter, the ages 
based on the Hipparcos parallax are slightly younger for 
t>8 Gyr. 

3.2. Sensitivity of the results 

Table 1 shows the results from four sets of calculations to 
check the sensitivity of the results caused by uncertainties 
in the input parameters. The mean age difference 5 gives 
the shift of the age due to variation of each of the input 
parameters. The variance a gives an indication about the 
scatter in the input parameter varied. 

All sets of calculations indicate that both 5 and a are 
significantly larger for the MS stars. This shows that the 
stars which contribute mostly to a large age spread are 
automatically identified as MS stars. The analysis further 
shows that the 15% uncertainty in the distance of the SGB 
stars is the main cause of their large age spread. An un- 
certainty of 5% is required to obtain variances comparable 
to those obtained for the other input parameters. A bet- 
ter determination of T cS and/or [Fe/H] does not lead to 
a significant improvement as long as the distances are not 
accurate down to a 5% level. 

Tabic 2 shows that a slightly larger improvement is ob- 
tained when the Hipparcos parallaxes are used. A reduc- 
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Fig. 7. (a) The AM R for the stars in Tables 7&8 with a 
dispersion in log(age) <0.05 (see Fig. 2 for the symbol descrip- 
tion) and (b) the AMR for the same stars, but now with the 
ages from Edv93ea. 



tion of the uncertainty in the ages will not be obtained 
from even more accurate distances, but from a better def- 
inition of the effective temperatures for the stars in the 
sample. Figure 7a further demonstrates that the number 
of stars with reliable ages has increased considerable with 
respect to the stars displayed in Fig. 4a. 

3.3. MS stars 

Figures 3 and 6 show the age-metallicity of all the stars 
in the sample analysed. The uncertainties in the ages of 
old MS stars are quite large. Near the MS the evolu- 
tion is quite slow and the isochrones of different ages are 
packed closely together. Small uncertainties in log T e g , My 
or metallicity give rise to large differences in the age of the 
MS star. In the age-metallicity plane this results in an in- 
creased scatter of the MS stars. 

Edv93ea removed from their sample stars which lied too 
close to the ZAMS with errors log(age) > 0.15. This is a 
subjective operation, because it partly depends on the ex- 
act definition of the ZAMS. As a consequence some MS 
might have been overlooked or too many taken out. 

The older ages obtained by Edv93ea for some of the 
solar metallicity MS stars (see Figs. 2 & 3) and the inclu- 
sion of these stars in the definition of the AMR results 
in a larger scatter in the definition of the disc AMR. The 
location of the MS turn-off tabulated by Bert94ea is used 
by us to distinguish a MS from a SGB star. In this way 
we classified more stars than Edv93ea as a MS star. The 



larger variance of the MS stars in Tables 1 & 2 appears to 
justify the current approach. 

As argued in Sect. 2.2.2 our fitting method is not substan- 
tially different from the one used by Edv93ea. The reason 
why they identified a different amount of stars too close 
to the MS is related to either of the following: 

— the exact location of the ZAMS adopted in the defini- 
tion of AMy for different metallicities, 

— a broader MS band due to convective overshoot. 

To some extent a different choice of the isochrone with 
metallicity Z from a specific [Fe/H] or [Mc/H] affects not 
only the age, but in our case also the distinction between 
MS or SGB star. 

3.4. SGB stars 

The situation changes if a star is on the SGB, where the 
evolution is relatively fast and the uncertainty in the age is 
considerably smaller. These stars potentially define a gen- 
uine disc AMR. Figure 4a shows the AMR for the SGB 
stars. It still has a large scatter in age and/or metallic- 
ity, but there is definitely a small slope of ^0.04 dex/Gyr 
present. For comparison we show in Fig. 4b&4c the AMR 
for the same stars, but now with respectively the visually 
assigned ages and the ages from Edv93ea. Figure 4b has 
a similar slope, while Fig. 4c shows an even steeper slope 
of ^0.07 dex/Gyr. The latter slope is also obtained when 
the Lutz-Kelker correction is applied to the absolute mag- 
nitude of the stars prior to the computation of the age or 
when the distances are derived from the Hipparcos par- 
allax, see Fig. 7a. However, one should be aware of the 
caveat that the slope might be shallower. Because some 
old, metal-rich stars might be absent in the sample, due 
to the selection criteria used by Edv93ea. 

Figures 3, 4a -c, 6 and 7a &b indicate that there is no 
apparent AMR for stars with an age less than 10 Gyr. 
We are basically dealing with a large metallicity spread 
among the stars. Only, if we consider the older stars a 
slope appears. 

A comparison of Figs. 4a&4b further indicates that 
there is no significant difference in the AMR slope be- 
tween the automatically determined ages and the visu- 
ally improved ages of the SGB stars. Preference is given 
to the objectively determined ages for the SGB stars in 
Fig. 4a, which have an uncertainty in their age of less than 
12%, while the estimated uncertainty in the metallicity is 
smaller than 0.07 dex. The slope in the AMR is not ex- 
pected to be an artifact. The large spread in age and/or 
metallicity in Figs. 4a and 7a&7b is likely real. It could 
imply that the AMR in Fig. 6 is a superposition of a mul- 
titude of relations due to in-fall or past mergers events. 

3.5. Initial stellar masses 

Figure || displays the ages with the corresponding masses 
obtained for the reliable (a T < 0.05) MS stars from Table 7 
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Fig. 8. The age-mass relation for the MS and SGB stars in 
Tables 7&8 with an uncertainty in log(age) < 0.05. The re- 
sults were obtained with calculations based on the Hipparcos 
parallax. We refer to Fig. 2 for the symbol description 

and the SGB stars from Table 8. One expects a smooth 
and continuous relation, but between 10-12 Gyr there is 
a hint of a slight discrepancy in the distribution of the 
stars. It could be an indication that stars with a larger age 
are younger and/or that the stars with ages larger than 
~ 8 Gyr are possibly slightly older. It might originate from 
the transition of isochrones with convective overshoot to 
isochrones without it, see Sect. 3.1 for the discussion about 
the presence or absence of overshoot. 
On the other hand, some stars with ages ^ 13 Gyr could 
be systematically younger. Since the same procedure was 
applied to determine the age, the origin should be due to 
one or more input parameters, possibly a combination of 
[Fe/H] and T e gf. However, [Fe/H] is likely the dominant 
factor. Moreover, at [Fe/H] £-0.4 Edv93ea found that 
the stars are relatively over-abundant in the a elements. 
In these cases we might have not used the optimum set of 
isochrones to determine the age, which results in a small 
discrepancy in both the older ages and the initial mass for 
these stars. 
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3.6. In summary 

We have computed with the Bert94ea isochrones new ages 
for the stellar sample defined by Edv93ea. The revised 
values are systematically larger towards older ages. The 
differences are considerably smaller when the distances 
are based on the Hipparcos parallax. The stars on the 
SGB define a disc AMR with a slope of ~ 0.07 dex/Gyr. 
A comparable slope is obtained, when the Edv93ea ages 
are adopted for the same stars. 
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Table 3. The age of F- and G-type on or very near the main sequence stars in the solar neighbourhood. The age is an average 
over the sensitivity analysis performed, see Sects. 2.2.3 & 3.2 for additional details. Column 1 gives the identifier of the star, 
columns 2&3 give the mean mass & the standard deviation for this mean value, columns 4&5 give the average age together 
with its standard deviation, and column 6 shows the number of good fits (out of the maximum 7, see Sect. 2.2.3) over which 
the average is obtained. In addition, n = \ indicates that 1 out of the 7 fits converged. The values in columns 2&4 refer to the 
good fit, while the standard deviation is in that case calculated over the 'good' (1) and 'second best' (6) estimates. A value 
n — indicates not a good fit was obtained in all 7 possible cases and the mean & standard deviation refer in this case to the 
'second best' value 
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Table 4. The age of F- and G-type sub-giant branch stars in the solar neighbourhood. See caption of Table 3 for a description 
of columns 1-6 and additional details. In addition, column 7 refers to the visually assigned ages 
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0.05 


9.440 


0.066 


7 


9.37 


HR 


4158 


1.08 


0.03 


9.793 


0.073 


6 


9.71 


HR 


4277 


1.06 


0.03 


9.816 


0.091 


6 


9.79 


HR 


4285 


1.06 


0.04 


9.869 


0.067 


7 


9.87 


HR 


4421 


1.23 


0.05 


9.540 


0.065 


6 


9.42 


HR 


4683 


1.10 


0.05 


9.718 


0.075 


7 


9.61 


HR 


4734 


1.11 


0.04 


9.822 


0.056 


7 


9.67 


HR 


4981 


1.38 


0.09 


9.426 


0.081 


7 


9.37 


HR 


4983 


1.11 


0.03 


9.562 


0.175 


5 


9.63 


HR 


4989 


1.14 


0.04 


9.623 


0.070 


3 


9.61 
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HR 
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0.067 


7 
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HR 
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HR 
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8fi 
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O 
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HR 
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1.21 


0.05 


9.450 


0.046 


4 


9.51 


HR 


TORS 

oy uo 


0.93 


0.03 


10.083 


0.042 


g 


10.10 


HR 
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1.26 
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u.uu 


q fi03 

y .uuo 
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U.U / o 


R 
u 


q 4q 
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HR 


6189 


1.04 


0.04 


9.829 


0.051 


g 


9.69 


HR 


6202 


1.38 


0.07 


9.406 


0.061 


5 


9.34 


HR 


6458 


0.88 


0.03 


10.169 


0.817 


I 


10.26 


HR 


6541 


1.29 


0.07 


9.529 


0.082 


7 


9.43 


HR 
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1.31 


0.03 
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7 
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7 
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HR 
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HR 
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9.512 
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4 
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HR 
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04 
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i 
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y .40 


HR 


7061 


1.37 


0.10 


9.449 
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7 


9.36 


HR 
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7 
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HR 
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04 
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7 
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HR 
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g 
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HR 
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g 
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HR 
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04 
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i 
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HR 


8027 


1.11 


0.04 


9.710 


0.062 


g 


9.66 


HR 


8041 


1.16 


0.05 


9.741 


0.061 


7 


9.60 


HR 


8077 

OU 1 1 


1.29 


1 

U. -LU 


q ^47 


0.124 


K 
O 


q 4q 

y .4y 


HR 


8354 


1.05 


0.04 


9.800 


0.056 


7 


9.66 


HR 


ouuo 


1.10 


0.01 


9.757 


0.031 


5 


9.66 


HR 
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1.21 
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u.uu 
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y .U-Lo 
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u.uyo 


u 
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y .4y 


HR 


8729 


1.10 


0.04 


9.847 


0.060 


7 


9.69 


HR 


880^ 
oouo 


1.36 


0.05 


9.329 


0.025 


7 


9.34 


I 1 IV 
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oooo 


1.21 


07 

U.U 1 


q TOO 

y .oyu 
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7 


q ^fi 

y.ou 


HR 
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J — LO 
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u.uu 


q fi14 
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U. -LUO 
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q ^7 
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HD 


9fi1 ^ 


1.11 


04 


q 730 
y. / ou 
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HD 


1 4q38 
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1 07 


03 
u.uo 


q 80Q 
y .ouy 


049 


4 


q 73 
y. i o 


HD 


18768 


0.85 


0.03 


10.191 


0.049 


6 


10.20 


HD 


38007 


0.92 


0.03 


10.108 


0.047 


6 


10.12 


HD 


68284 


0.94 


0.04 


10.016 


0.069 


7 


10.02 


HD 


78558 


0.86 


0.02 


10.180 


0.034 


3 


10.20 


HD 


130551 


0.97 


0.02 


9.918 


0.036 


7 


9.82 


HD 


144172 


1.12 


0.04 


9.723 


0.048 


7 


9.61 


HD 


157089 


0.88 


0.03 


9.498 


0.564 


1 


10.25 


HD 


159307 


1.07 


0.04 


9.775 


0.067 


7 


9.78 


HD 


188815 


0.94 


0.03 


9.930 


0.072 


7 


9.93 


HD 


198044 


1.03 


0.04 


9.865 


0.064 


6 


9.82 


HD 


200973 


1.17 


0.05 


9.650 


0.065 


7 


9.51 


HD 


201099 


0.90 


0.02 


10.100 


0.035 


6 


10.14 


HD 


205294 


1.13 


0.03 


9.726 


0.043 


5 


9.58 


HD 


221830 


0.86 


0.03 


10.195 


0.046 


6 


10.20 
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Table 5. Similar to Table 3, but now with the Lutz-Kelker (1973) correction (n = 4, oy/7r = 0.15; Hanson 1979, equation 
applied to the absolute magnitude of the F- and G-type stars on or near the main sequence prior to the calculation of the 



ID 


m 


0~m 


r 


<7 T 


n 


HR 


140 


1.36 


0.05 


9.344 


0.030 


7 


HR 


219 


0.96 


0.03 


9.814 


0.197 


6 


HR 


370 


1.27 


0.05 


9.495 


0.037 


7 


HR 


458 


1.35 


0.05 


9.419 


0.030 


7 


HR 


672 


1.34 


0.06 


9.487 


0.051 


5 


HR 


784 


1.23 


0.03 


9.298 


0.289 


7 


HR 


962 


1.39 


0.05 


9.433 


0.033 


7 


HR 


1010 


0.96 


0.03 


9.952 


0.091 


7 


HR 


1173 


1.50 


0.05 


9.158 


0.046 


7 


HR 


1257 


1.50 


0.06 


9.334 


0.051 


7 


HR 


1489 


1.28 


0.07 


9.542 


0.088 


6 


HR 


1687 


1.48 


0.05 


9.115 


0.073 


7 


HR 


2220 


1.38 


0.04 


9.215 


0.087 


7 


HR 


2721 


0.98 


0.03 


9.776 


0.304 


7 


HR 


2943 


1.59 


0.06 


9.207 


0.031 


7 


HR 


3018 


0.80 


0.04 


8.754 


0.285 





HR 


3578 


0.80 


0.02 


10.176 


0.020 


2 


HR 


3954 


1.49 


0.06 


9.314 


0.033 


7 


HR 


4012 


1.46 


0.06 


9.379 


0.035 


7 


HR 


4067 


1.48 


0.06 


9.305 


0.031 


7 


HR 


4395 


1.59 


0.06 


9.230 


0.036 


7 


HR 


4529 


1.44 


0.05 


9.395 


0.035 


7 


HR 


4533 


1.42 


0.05 


9.258 


0.042 


7 


HR 


4540 


1.37 


0.05 


9.410 


0.029 


7 


HR 


4688 


1.44 


0.05 


9.296 


0.030 


7 


HR 


4785 


0.97 


0.03 


9.916 


0.136 


7 


HR 


4845 


0.81 


0.03 


10.242 


0.017 


2 


HR 


4981 


1.50 


0.06 


9.334 


0.044 


7 


HR 


5011 


1.29 


0.07 


9.522 


0.075 


5 


HR 


5235 


1.56 


0.08 


9.316 
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3 


HR 


5323 


1.43 


0.05 


9.379 


0.031 
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9 
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7 


HR 
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9 
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7 


HR 
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1.55 


0.10 


9 
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7 


HR 


6243 


1.73 


0.05 


9 
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0.034 


6 


HR 


6409 


1.60 


0.06 


9 


223 


0.037 


7 


HR 


6907 


1.41 


0.05 


9 


318 


0.028 


7 


HR 


7061 


1.51 


0.06 


9 


327 


0.037 


7 


HR 


7126 


1.48 


0.04 


9 


074 


0.130 


7 


HR 


7560 


1.34 


0.05 


9 


445 


0.034 


7 


HR 


7955 


1.40 


0.05 


9 


396 


0.030 


7 


HR 


8181 


0.90 


0.01 


9 


990 


0.051 


3 


HR 


8472 


1.61 


0.05 


9 


268 


0.036 


6 


HR 


8885 


1.45 


0.05 


9 


274 


0.027 


7 


HD 


6434 


0.86 


0.02 


10 
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0.054 


4 


HD 
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10 
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z 


HD 
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0.02 


10 


155 
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HD 
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0.79 


0.02 


10 


246 


0.000 


2 


HD 


78747 


0.80 


0.02 


10 


192 


0.000 


2 


HD 


98553 


0.90 


0.03 


10 


008 


0.158 


7 


HD 


114762 


0.82 


0.03 


10 


203 


0.831 


1 


HD 


134169 


0.82 


0.03 


10 


209 


0.056 


2 


HD 


148816 


0.84 


0.02 


9 


728 


0.652 


1 


HD 


174912 


0.86 


0.01 


10 


123 


0.044 


3 


HD 


199289 


0.82 


0.02 


9 


866 


0.039 


1 


HD 


201891 


0.79 


0.04 


9 


856 


0.364 


3 


HD 


208906 


0.83 


0.03 


10 


082 


0.140 


6 


HD 


210752 


0.84 


0.02 


10 


158 


0.023 


3 


HD 


215257 


0.87 


0.02 


10 


061 


0.011 


2 



Table 6. Similar to Table 5, but now for F- and G-type stars on the sub-giant branch 



ID 


m 


Cm 


r 


ay 


n 


HR 


17 


1.10 


0.04 


9.749 


0.061 


5 


HR 


33 


1.14 


0.06 


9.724 


0.072 


4 


HR 


35 


1.33 


0.04 


9.322 


0.039 


6 


HR 


107 


1.26 


0.06 


9.516 


0.093 


6 


HR 


145 


1.18 


0.06 


9.677 


0.092 


6 


HR 


203 


1.05 


0.04 


9.906 


0.057 


6 


HR 


235 


1.15 


0.02 


9.577 


0.116 


6 


HR 


244 


1.35 


0.06 


9.491 


0.060 


7 


HR 
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0.07 
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3 


HR 
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1.30 


0.07 
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9.342 
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HR 
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9.556 


0.063 


7 


HR 


483 


1.08 


0.03 


9.842 


0.054 


7 


HR 


573 


1.10 


0.04 


9.736 


0.065 


6 


HR 


646 


1.29 


0.05 


9.501 


0.046 


4 


HR 


720 


1.06 


0.04 


9.884 


0.064 


7 


HR 


740 


1.45 


0.10 


9.366 


0.087 


7 


HR 


799 


1.23 


0.03 


9.369 


0.164 


7 


HR 


1083 


1.48 


0.05 


9.251 


0.022 


7 


HR 


1101 


1.14 


0.04 


9.761 


0.054 


4 
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Table 6. Continued ... 



ID 


m 


0~m 


T 


ay 


n 


HR 


1294 


1 


03 


0.05 


9.950 


0.069 


6 


HR 


1536 


1 


34 


0.06 


9.532 


0.062 


6 


HR 


1545 


1 


21 


0.04 


9.583 


0.076 


6 


HR 


1673 


1 


31 


0.08 


9.471 


0.091 


7 


HR 


1729 


1 


14 


0.03 


9.784 


0.046 


6 


HR 


1780 


1 


18 


0.03 


9.504 


0.127 


7 


HR 


1983 


1 


25 


0.03 


9.357 


0.141 


7 


HR 


2047 


1 


08 


0.03 


9.421 


0.662 


5 


HR 


2141 


1 


01 


0.03 


9.944 


0.040 


6 


HR 


2233 


1 


32 


0.07 


9.487 


0.072 


7 


HR 


2354 


1 


17 


0.05 


9.729 


0.063 


7 


HR 


2493 





98 


0.03 


9.932 


0.054 


5 


HR 


2530 


1 


24 


0.06 


9.522 


0.087 


6 


HR 


2548 


1 


35 


0.09 


9.437 


0.086 


7 


HR 


2601 


1 


05 


0.04 


9.838 


0.068 


7 


HR 


2835 





97 


0.03 


9.921 


0.046 


6 


HR 


2883 





89 


0.02 


10.078 


0.045 


6 


HR 


2906 


1 


32 


0.06 


9.520 


0.059 


7 


HR 


3176 


1 


23 


0.06 


9.673 


0.075 


7 


HR 


3220 


1 


32 


0.05 


9.428 


0.067 


7 


HR 


3262 


1 


15 


0.03 


9.548 


0.120 


3 


HR 


3271 


1 


31 


0.06 


9.557 


0.064 


5 


HR 


3538 


1 


03 


0.03 


9.946 


0.044 


7 


HR 


3648 


1 


15 


0.03 


9.772 


0.048 


6 


HR 


3775 


1 


33 


0.10 


9.483 


0.092 


6 


HR 


3881 


1 


19 


0.04 


9.706 


0.055 


5 


HR 


3951 


1 


26 


0.05 


9.631 


0.064 


6 


HR 


4027 


1 


16 


0.06 


9.758 


0.079 


6 


HR 


4039 


1 


03 


0.04 


9.855 


0.072 


6 


HR 


4150 


1 


42 


0.08 


9.376 


0.071 


7 


HR 


4158 


1 


11 


0.03 


9.761 


0.038 


7 


HR 


4277 


1 


08 


0.04 


9.833 


0.082 


7 


HR 


4285 


1 


11 


0.04 


9.796 


0.056 


6 


HR 


4421 


1 


29 


0.06 


9.472 


0.074 


6 


HR 


4657 





92 


0.02 


9.926 


0.079 


7 


HR 


4683 


1 


14 


0.05 


9.676 


0.059 


6 


HR 


4734 


1 


17 


0.05 


9.742 


0.067 


7 


HR 


4767 


1 


09 


0.02 


9.689 


0.110 


6 


HR 


4903 


1 


46 


0.06 


9.416 


0.068 


6 


HR 


4983 


1 


14 


0.03 


9.600 


0.072 


6 


HR 


4989 


1 


17 


0.04 


9.643 


0.074 


7 


HR 


5019 


1 


00 


0.03 


10.028 


0.048 


6 


HR 


5338 


1 


44 


0.06 


9.420 


0.050 


5 


HR 


5353 


1 


43 


0.12 


9.468 


0.102 


3 


HR 


5423 


1 


33 


0.01 


9.557 


0.009 


4 


HR 


5447 


1 


30 


0.04 


9.357 


0.033 


7 


HR 


5459 


1 


04 


0.03 


9.912 


0.046 


5 


HR 


5691 


1 


31 


0.06 


9.524 


0.065 


7 


HR 


5868 


1 


13 


0.04 


9.780 


0.056 


7 


HR 


5914 





89 


0.03 


10.115 


0.044 


6 


HR 


5933 


1 


24 


0.06 


9.528 


0.098 


6 


HR 


5968 





97 


0.04 


10.037 


0.055 


7 


HR 


5996 


1 


35 


0.08 


9.521 


0.073 


6 


HR 


6189 


1 


11 


0.05 


9.744 


0.065 


7 



ID 


m 


0~m 


T 


a T 


n 


HR 


6202 


1.45 


0.09 


9.354 


0.074 


5 


HR 


6315 


1.16 


0.03 


9.512 


0.136 


7 


HR 


6458 


0.89 


0.03 


10.162 


0.047 


6 


HR 


6541 


1.33 


0.05 


9.491 


0.049 


5 


HR 


6569 


1.38 


0.04 


9.347 


0.045 


7 


HR 


6598 


1.03 


0.05 


9.918 


0.080 


7 


HR 


6649 


1.04 


0.04 


9.864 


0.046 


7 


HR 


6701 


1.35 


0.05 


9.482 


0.051 


5 


HR 


6775 


0.90 


0.01 


10.030 


0.045 


2 


HR 


6850 


1.32 


0.07 


9.410 


0.083 


7 


HR 


7232 


1.04 


0.04 


9.957 


0.063 


7 


HR 


7322 


1.20 


0.06 


9.620 


0.095 


6 


HR 


7534 


1.31 


0.05 


9.451 


0.059 


7 


HR 


7766 


1.00 


0.04 


9.940 


0.052 


7 


HR 


7875 


1.06 


0.04 


9.847 


0.052 


6 


HR 


8027 


1.13 


0.03 


9.711 


0.044 


7 


HR 


8041 


1.23 


0.05 


9.656 


0.067 


7 


HR 


8077 


1.32 


0.09 


9.529 


0.098 


7 


HR 


8354 


1.11 


0.05 


9.720 


0.069 


6 


HR 


8665 


1.16 


0.05 


9.693 


0.080 


6 


HR 


8697 


1.26 


0.05 


9.566 


0.065 


4 


HR 


8729 


1.17 


0.05 


9.757 


0.071 


7 


HR 


8805 


1.42 


0.07 


9.340 


0.069 


7 


HR 


8853 


1.24 


0.06 


9.585 


0.092 


7 


HR 


8969 


1.23 


0.04 


9.583 


0.073 


5 


HD 


2615 


1.17 


0.05 


9.640 


0.073 


7 


HD 


14938 


1.10 


0.04 


9.771 


0.052 


7 


HD 


17548 


0.87 


0.02 


10.111 


0.040 


7 


HD 


18768 


0.89 


0.04 


10.115 


0.074 


7 


HD 


30649 


0.86 
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Table 7. The age of F- and G-type on or very near the main sequence stars in the solar neighbourhood. The age is an average 
over the sensitivity analysis performed, see Sects. 2.2.3 & 3.2 for additional details. The distances of the stars are based on the 
Hipparcos parallax (ESA 1997). Column 1 gives the identifier of the star, columns 2&3 give the mean mass & the standard 
deviation for this mean value, columns 4&5 give the average age together with its standard deviation, and column 6 shows the 
number of good fits (out of the maximum 7, see Sect. 2.2.3) over which the average is obtained. In addition, n = 1 indicates 
that 1 out of the 7 fits converged. The values in columns 2&4 refer to the good fit, while the standard deviation is in that 
case calculated over the 'good' (1) and 'second best' (6) estimates. A value n — indicates not a good fit was obtained in all 7 
possible cases and the mean & standard deviation refer in this case to the 'second best' value 
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Table 8. The age of F- and G-type sub-giant branch stars in the solar neighbourhood. See caption of Table 7 for additional 
details 
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